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Abstract. We construct low regularity solutions of the vacuum Einstein constraint equations on compact 
manifolds. On 3-manifolds we obtain solutions with metrics in H s where S > 3/2. The constant mean 
curvature (CMC) conformal method leads to a construction of all CMC initial data with this level of regularity. 
These results extend a construction from [Ma04] that treated the asymptotically Euclidean case. 



1. Introduction 

Initial data for the vacuum Cauchy problem in general relativity consists of a Riemannian 
manifold (M 3 , g) and a symmetric (0, 2)-tensor K. One seeks an isometric embedding of 



(M, g) into a Ricci-flat Lorenzian manifold such that K is the extrinsic curvature of this 



embedding. Since the ambient manifold is Ricci flat, the initial data (M, g, K) must satisfy 
the compatibility conditions 

7" R g - \K\l + ti g K 2 = 

^ : 9 (!) 

bjr). div 5 K - d tig K = 0, 

where R g is the scalar curvature of the metric g. System (1) is known known as the Einstein 
constraint equations; the first equation is the Hamiltonian constraint and the second is the 
momentum constraint. 

The Cauchy problem has a formulation as a second-order quasilinear system of hyperbolic 
PDEs called the Einstein equations. In [CB52] Choquet-Bruhat proved that this problem 
is locally well-posed for sufficiently smooth initial data satisfying the constraint equations. 
Subsequent investigations have reduced the amount of regularity required from the initial 
data, but development of the low regularity theory of the constraint equations has lagged 
behind that of the evolution problem. For example, the local well-posedness theorem of 
[HKM77] for the Einstein equations permits initial data (g, K) in Hf x Hf -1 with s > 5/2. 
But solutions of the constraint equations with this level of regularity were not known to 
exist until recently [CBIY00]. The need for a theory of rough solutions of the constraint 
equations became especially prominent after the appearance of [KR] which proved an a 
priori estimate for the the time of existence of solutions of the vacuum Einstein equations 
depending on the size of (dg, K) in _£P -1 x H 13 ^ 1 with s > 2. 
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Constructions of rough solutions of the constraint equations have been obtained in [CB04] 
on compact manifolds and in [CBIYOO] and more generally in [Ma03] on asymptotically 
Euclidean manifolds. These works provide solutions with metrics in H 2 . More specifically, 
these solutions were constructed in W^' p spaces where k > 2 is an integer and k > 3/p. 
The restriction k > 3/p is interesting because it is analogous to condition s > 3/2 for 
H s spaces. Since s = 3/2 is the scaling limit for the Einstein equations, this level of 
regularity is a natural lower bound for for seeking solutions of the constraint equations. 
These considerations motivated work in [Ma04] which gave a construction of all constant 
mean curvature (CMC) asymptotically Euclidean solutions of the constraint equations with 
metrics in Hf for real values of s > 3/2. In this paper we extend the low regularity 
construction from [Ma04] to compact manifolds. In particular, we show that the conformal 
method can be used to construct all CMC vacuum solutions (g, K) of the constraint equations 
belonging to H s x H s ~ l with s > 3/2. 

The constant mean curvature conformal method of Lichnerowicz [Li44] and Choquet- 
Bruhat and York [CBY80] is the basis of all known constructions of rough solutions of the 
constraint equations. The method starts with conformal data (M ra , g, S, r) where M n is a 
smooth, compact, connected manifold with n > 3, g is a metric on M, S is a trace-free 
(0, 2)-tensor, and r is a scalar function. 1 We then seek a solution (g, K) of the constraint 
equations on M of the form 



g = <t> 2K g 



K = <p~ 2 (a + LW) + -g 

n 

where the unknowns (0, W) are a positive function and vector field respectively, k = 
is a dimensional constant, and L is the conformal Killing operator defined by 

LW ab = V a W h + V h W a - - div Wg ab . 

n 



It follows that (g, K) solves the constraint equations so long as 

2 rh 2K+1 

(2) 



-aA g( j) + R g( p = \S + LW\ 2 g (f)~ 2K - 3 - br 2 2K+1 



divW = 6</r nK dr-divS 

where R g is the scalar curvature of g and a = 2k + 4 and b = axe constants. The 
function r specifies the mean curvature tr^ K, and the conformal method is most successful 
when t is assumed to be constant (see, however, [MI96][CB04][IOM04] for results in 
the non-CMC case). When r is constant, system (2) decouples and the second equation 
becomes 

divhW = -divS (3) 



x We work with general n-dimcnsional manifolds since some of our auxiliary results are meaningful outside the context 
of the constraint equations. 
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In the smooth setting, equation (3) is always solvable and leads to a decomposition of the 
set of smooth symmetric, trace-free (0, 2)-tensors 

y = imL©kerdiv (4) 

known as York splitting. Tensors in D ker div are known as transverse-traceless, and 
specifying conformal data (M,g,S,r) for the CMC conformal method is equivalent to 
specifying (M, g, a, r), where a is the transverse-traceless part of S. The study of system 
(2) reduces to determining conditions when the Lichnerowicz equation 

-aA g <j) + R g <f) = \a\ 2 g 0- 2k ~ 3 - 6t 2 2k+1 (5) 

has a positive solution. For smooth metrics it is known[CBY80][Is95] that (5) is solvable 
if and only if one of the following conditions hold: 

1. cr^Oandr^O, 

2. X g > 0, a ^ and r = 0, 

3. X g = 0, o = andr = 0, 

4. Xg < 0, o = and r ^ 0. 
Here X g is the Yamabe invariant defined by 

. , j M *\Vf\l + R 9 f 2 dVg 

X n = mf 



9 /GC°°(M) " ~ II/M2 



f& UJ L 2 * 

Our main result, Theorem 4.5, shows that these conditions are also necessary and sufficient 
for the solvability of (5) when g e H s and a e H s ~ l with s > n/2. We also show in 
Theorem 5.2 that the York decomposition (4) holds for symmetric, trace-free (0, 2)-tensors 
in H 3 ^ 1 and hence the CMC conformal method fully extends to the rough setting. 

Aside from lowering the regularity of solutions, we also obtain some technical improvements 
to the previous construction of rough solutions on compact manifolds obtained in [CB04]. 
Solutions were found in that reference in the Yamabe negative case (X g < 0) only under 
an additional hypothesis on the singularities of a, namely a e -fP -1 H L°°. Also, the case 
X g > 0, r = and a ^ when a has zeros was not treated. The zeros and singularities 
of a pose difficulties for the method of sub- and supersolutions used to find solutions 
of (5). Our approach to finding sub- and supersolutions of (5), discussed in Section 4, 
circumvents these problems by means of an initial conformal transformation to a metric 
with the singularities and zeros of a already encoded in the scalar curvature. This strategy 
has been used previously for smooth conformal data, but the choice of conformal change 
used here is new and has interest even in the context of smooth solutions. 



4 



David Maxwell 



2. Considerations for Rough Metrics 

Before turning to the specifics of solving the constraint equations, we start with an outline of 
the rules for working with rough metrics. Throughout this paper, let (M n , g) be a compact, 
connected Riemannian manifold with n > 3 and g E H s . Our definitions and notation for 
Sobolev spaces follows [Ta96]. The value s — n/2 is a dividing point between metrics that 
are natural for geometric analysis and those that are not. For example, the Laplacian of a 
metric in H s with s < n/2 doesn't necessarily take H 2 to L 2 or even H 1 to if -1 . This is a 
consequence of the fact that a metric in H s with s < n/2 is not necessarily in L°°. On the 
other hand, a metric in H s with s > n/2 can be nearly as good as a smooth metric. The 
following well-known lemma is of central importance. 

Lemma 2.1. Suppose M n is a compact manifold, a < min(si,S2), s\ + S2 > and 
er < s\ + S2 — § • Then pointwise multiplication extends to a continuous bilinear map 



Lemma 2.1 implies that if u E H s with s > n/2, then multiplication by u takes H a to 
H a for every a E [—s, s]. In particular, H s is an algebra. This fact allows us to have well 
defined conformal changes of metric, for example. If is a positive function in H s , and if g 
is a metric in H s , then Lemma 2. 1 implies 4> a g is again rough metric in H s for any positive 
integer a. To allow arbitrary real exponents a, and to work with other nonlinearities, we 
require the following consequence of Lemma 2.1. 

Lemma 2.2. Suppose M n is a compact manifold and F : I ^ R. is smooth on some interval 
I. Ifs > n/2 and u E H S (M) satisfies u(x) E I for all x E M, then F(u) E H S (M). 

Proof: Working in a coordinate chart we can assume that u is defined on R n and need only 
show that F(u) E Hf . Let k be the largest integer such that k < s. If a is a multi-index 
such that |a| = k, then d a F{u) is a sum of terms of the form F^ r \u)d ai u ■ ■ ■ d ar u where 
Jj\ = \ \oti\ = k. From repeated applications of Lemma 2.1 we have d ai u ■ ■ ■ d ar u E H° qc 
for any a < s - k + (r - l)(s - f ). Since F r {u) E this shows that F(u) E if£ c . To 
obtain the full regularity of F(u), consider any derivative diF(u) = F'(u)diU. Since F' is 
smooth, we have F'(u) E H^ qc and from Lemma 2.1 we have F'{u)diU E H ai ~ 1 for every 
o\ < s such that o\ < k + s — ^. Hence F(u) E H^ c - Repeating this argument shows 
that F(u) E for every 02 < s such that 02 < + 2 (s — |). A finite number of further 
repetitions shows that F(u) E Hf . □ 

Consider the Laplacian A g of a rough metric g. In local coordinates we have 



H Sl (M) x i/ S2 (M) -> H a (M). 
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It follows from Lemmas 2.1 and 2.2 that the Laplacian has the form in local coordinates 

a ij didj + h'i); 

where a e and 6* e if* -1 . Using Lemma 2.1 it follows that A p : — > iJ " -2 for 
any a G [2 — s, s]. The following class of differential operators encodes the regularity of 
the coefficients of the Laplacian. 

Definition 2.3 Let A be a linear differential operator that in any coordinate chart has the 
form 

A= 2 a a d a , 

\a\<m 

where a a is aR kxk valued function. We say that 

A G £ ?m > s 

ifa a G H°~ m+lal for all \a\ < m. 

An analogous definition holds for differential operators on sections of vector bundles. From 
Lemma 2.1 we have the following mapping properties. 

Lemma 2.4. Suppose m < 2s and o G [m — s, s]. If A e J? rn,s , then A defines a 
continuous map 

A : H a — > iJ a - m . 

We say A is elliptic if for each x, the constant coefficient operator 2 | a |=m flQ (^^a is elliptic 
in the usual sense. Using Lemmas 2.1 and 2.2 it is easy to check that the Laplacian A g , 
conformal Laplacian —aA g + R g , and vector Laplacian Al = divL all belong to Jzf 2 ' s 
when s > n/2 and are all elliptic. 

The following interior elliptic estimate was proved in [Ma04]. We use the notation A < B 
to mean A < cB for a positive constant c where c does not depend on the variable functions 
and parameters of A and B. For example, in estimate (6) below the implicit constant is 
independent of the function u but can can depend on the open sets U, V, the operator A, 
and so forth. 

Proposition 2.5. Let U and V be open sets in R n with U CC V, and suppose s > n/2 and 
a E (m — s, s] . If A e m,s is elliptic, then for every u e H a we have 

\\u\\h°(U) Z \ \Au\\ H <r-m(y) + \ H\H m - s {V)- (6) 

Using coordinate charts and cutoff functions, we obtain as an immediate consequence 
a-priori estimates on compact manifolds. 
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Corollary 2.6. Suppose (M n ,g) is a compact Riemannian manifold with g E H s and 
s > n/2. If A is a second order elliptic differential operator in J2? 2 ' s , then 

\\ x \\h s (M) ^ \\AX\\ H 8-2( M ) + ||X|| L 2( M ) 

for all X e H S (M). 

We now turn to duality properties and integration by parts for rough metrics. We would 
like to claim, for example, that the equation 

I (— aA a u + R a u)v dV a = \ a \Vu\i + R a uv dV a 
JM JM 9 

that holds for the conformal Laplacian of a smooth metric also holds for rough metrics. We 
need to show that integration by parts holds and that the term J M R g uv dV g can be suitably 
interpreted for rough metrics (note that R g need not be a locally integrable function when 

n = 3). 

Suppose g is a rough metric in H s with s > n/2 and g is a smooth background metric on 
M. We have a bilinear form defined for u,v £ C°°(M) by 

( u i v )(M,g) = f M uvdV 9 ( 7 ) 
which extends for any a e R to a continuous bilinear map 

(;-) {M jy.H°xH-°(M)^R 

and induces an isomorphism 

H a (M) = (H- a (M))* . 

In applications we would prefer use the rough metric g instead of the smooth metric g to 
induce this isomorphism. This can be done for a restricted range of Sobolev spaces. From 
Lemma 2.2, we have dV g = h dV g where h > and h E H S (M). Multiplication by h is 
continuous on on H a for every a G [— s, s], so we can define for u E H a and v E H~ a 

( U ' V )(M,g) = ( hu , v )(M,g)- 

With this definition, we have the following easy consequence. 

Lemma 2.7. Suppose (M n , g) is a compact Riemannian manifold with g E H s ands > n/2 
and suppose a E [—s, s\. Then for every u E H a and v E H~° , 

\( U ' V )(M,g)\ ^ \\u\M\vWh-'- 
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Moreover, H a (M) is naturally isomorphic to CT (M))* through the pairing (•, -)(Mg)- 
Finally, for smooth functions u and v, 

( U > V )(Mj)=f M ™ dV 9- 

These ideas readily extend to sections of vector bundles as well. We show here explicitly 
how this is done for sections of the tangent bundle. The smooth background metric g gives 
rise to a bilinear form defined for smooth vector fields X and Y by 

This map extends to a continuous bilinear form of less regular vector fields 

for every a G R. Now if X is a H a vector field with a G [— s, s], we can define X* G i/' 7 
by 

where as before dV g = h dV g . We then define for vector fields X and Y in H a and iJ _cr 
respectively 

and we note that if X and Y are smooth then 

To establish an integration by parts formula, suppose that u G H a with a G [2 — s, s], 
f G i/ 2_cr , and suppose that v is supported in the domain of some coordinate chart. We 
then have in local coordinates 

(Vu, Vw) (M;5) = J] Udetgg lJ diU, djv) 

i,j ' '(R-™)<?) 

where the bilinear form on the right hand side is the usual pairing of against i/ 1_cr 

for the flat Euclidean metric. The standard integration by parts formula 

is valid for / G H T (W l ) and g compactly supported in iJ 1_r (R. n ), for any r G R. Hence 
2 Udet gg lJ diu, djv) = J] (-djJdet gg tJ diU, v) = (-A g u, v), M y 

A partition of unity argument then establishes 

(Vu,Vu) (M>fl) = (-A s u,u) (M>9) 

for every m G H a (M) and t> G H 2 ~ a (M) with cr G [2 — s, s]. A similar approach also 
shows that integration by parts also holds for the vector Laplacian. 
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Lemma 2.8. Suppose (M n ,g) is a compact Riemannian manifold with g G H s with 
s > n/2. Then for functions u G H a and v G i/ 2 ~°" with a G [2 — s, s] we have 

(Vu,Vv) {Mjg) = (-A g u,v) {Mgy 

For vector fields X G H a andY G H 2 ~ a with a G [2 - s, s] we have 

(LX,Lr) (Mi3) = (-A L X,F) (M)9) . 

Finally, we show that maximum principles hold for rough metrics. These are used most 
importantly to ensure that the conformal factors we construct are positive. We say that a 
distribution V G H a with a G [— s, s] satisfies V > if 

(va)(m, 9) >o 

for every smooth function <\> > 0. 

Lemma 2.9. Suppose (M n , g) is a compact Riemannian manifold with g G H s and 
s > n/2. Suppose also that V G H s ~ 2 and that V > with V ^ 0. If u G H s satisfies 

-Au + Vu < 0, (8) 

then u < 0. 

Proof: Let v = . One readily finds from Lemma 2.1 , since G H 1 , that uv belongs 
to H 1 . Since s — 2 > —1 and since V G H s ~ 2 , we can apply V 7 to Since V > 0, 

and since itt) > 0, we have (V, uv)/ Mg \ > 0. Finally, since u satisfies (8) we obtain 

(-Au,v) {Mjg) < - (V,uv) {Mjg) < 0. 

Integrating by parts and using the fact that 

Vi> = { x y = 1 

we conclude v is constant. To show that v = we see that if not, then = v is a 
positive constant and hence a is a positive constant. Since —Am + Vu = Vu > 0, and 
since —Au + Vu < Owe conclude V < and V > 0. But then V = 0, a contradiction. 

□ 

A version of the following strong maximum principle was proved in [Ma04] for asymp- 
totically Euclidean metrics. Since the argument there was purely local, using only the 
regularity of the metric and the connectivity of the manifold, the same proof yields a strong 
maximum principle for compact manifolds. 

Proposition 2.10. Suppose (M n , g) is a compact Riemannian manifold with g G H s and 
s > n/2. Suppose that u G H S (M) is nonnegative and that V G H S ~ 2 (M) satisfies 

—A g u + Vu > 0. 

Ifu(x) = at some point x G M, then u vanishes identically. 
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3. The Yamabe Invariant of a Rough Metric 

The Yamabe invariant for smooth compact manifolds, defined by 



J M a\Vf\ 2 g + R g p dV g 



\,, = inf — 

9 /ec°°(M) II/M2 



where 2* = 2n/(n— 2), plays an essential role in classifying CMC solutions of the constraint 
equations on compact manifolds. To understand properties of the Yamabe invariant for 
rough metrics, we consider the quadratic bilinear form 



A(f, h)=f M a (Vf, Vh) + R g fh dV g . 



When g is smooth then A is obviously continuous on H 1 x H 1 . This is also true when 
g is rough, so long as we interpret J M R g fh dV g as (R g , fh), Mg y To see this, we note 
that the term J M a (V/, Vh) dV g poses no difficulties and we need only consider the scalar 
curvature term. From Lemmas 2.1 and 2.8 we have 

(R 9 Jh} {M J<\\R 9 \\ H ^\\f\\m\\h\\m (9) 

since 1 + 1 — ^ > 2 — s. This establishes the continuity of A on H 1 . In fact, we can 
say more than this. Pick 77 such that < 77 < 1 and such that 2r? < s — f • Then 
(1 — rj) + (1 — rj) — > 2 — s as well, and we conclude from Lemma 2.1 that 

( R g,f h )(M, 9 ) ^\\ R g\\H^\\f\\H^-n\\h\\ H i-v (10) 
Hence the scalar curvature term is continuous on H 1 ^ 11 x H 1 ^ 1 '. 

We would like to show that \ g ^ —00, which we do by considering the first eigenvalue of 
the conformal Laplacian. Let 



where 1 1/| \u> is the metric dependent norm defined by 



L = 1,1/1" "v, 



We have 



\„ = inf J 2 * (/) 



and the first eigenvalue of the conformal Laplacian is defined by 

Uh = inf hif)- 

fro 

If Xg < then it is easy to see that ji g < X g and therefore X g 7^ —00 if ji g 7^ —00. Our 
estimate for a lower bound for ii g starts with the following estimate for A. 
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Lemma 3.1. Let (M n , g) he a compact Riemannian manifold with g E H s and s > n/2. 
Then there are positive constants C\ and C2 such that for every f E H 1 

A(f,f)>C 1 \\f\\ 2 Hl -C 2 \\f\\l> 

Proof: Pick rj as was done in deriving inequality (10). From interpolation we have for any 
e > that 

\\f\\m-,<C(v,e)\\f\\ L i + e\\f\\ m . 
From inequality (10) we then have 



We also note that 



Since 



+ 



2 L2 + e 2 \\f\\ 2 H , 



L 2 ' 



(11) 



(12) 



A(fJ)> J M a\Vf\ 2 g dV g -\(R g ,f\ Mg) 
we obtain from (11) and (12), taking e sufficiently small, that 

A(f,f) > \\f\\ 2 Hl - (C(r),e,\\R g \\ HS -2) + 1) ||/||| 2 . 

□ 

It follows immediately that X g and \i g are both finite. 

Lemma 3.2. Let (M n , g) he a compact Riemannian manifold with g E H s and s > n/2. 
Then X g and fi g are both finite. 

Proof: Since ||/||^ 2 ^ ||/||^ 2 *, it follows that if A 3 = — 00 then \i g = —00. So itis enough 
to show that fi g is finite. From Lemma 3.1 we have 

j M a\Vf\ 2 g + R g f 2 ^^Cill/H^-Csll/lg^-^H/lg,. 
Dividing by 1 1/| |^ 2 we conclude that 

/A? > -C 2 . 



□ 
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Proposition 3.3. Let (M n , g) he a compact Riemannian manifold with g e H s and s > 
n/2. Then the following are equivalent: 

1 . There exists a positive function <fi e H s such that g = § 2n g satisfies Rg is 
continuous and positive (resp. zero, resp. negative). 

2. \ g is positive (resp. zero, resp. negative). 

3. /jg is positive (resp. zero, resp. negative). 

Proof: We start by proving Condition 3 implies Condition 1 . We claim that there is function 
/ in H 1 , not identically zero, such that foif) = l^g- To see this, let {/&} be a sequence of 
smooth functions such that | \f k \ \ H i = 1 and hifk) converges to /j, g . Since the sequence is 
bounded in H 1 , we can assume without loss of generality that it converges weakly to some 
/ G H 1 . From Lemma 3.1 we have 

r/f w C lWfk\\ 2 H i -C 2 \\fk\\ 2 L 2 Ci 

Uh) - — \m — = 

Since — > / strongly in L 2 , and since the sequence J2(fk) is bounded, we conclude that / 
is not identically zero. It is easy to see that J2 is weakly lower semicontinuous on H 1 \ {0} 
and hence we conclude that J2U) = /V 

Since / is a minimizer in H 1 of J 2 (/) it satisfies 

-ttAgf + Rgf = flgf 

in H^ 1 . Proposition 6.1 then implies / e H s . We can assume without loss of generality 
that / > since |/| e H 1 and J2(f) = -hdfl)- So from the strong maximum principle 
Proposition 2.10 we conclude that / > everywhere. Letting 4> = f and g = (p 2K g, it 
follows that Rg = _2k_1 (— a A g (p + R g 4>) = <\T 2K fig Since <p is continuous and positive, 
we conclude that Rg is continuous and everywhere has the same sign as pig. 

We next show Condition 1 implies Condition 2. Since \ g is a conformal invariant, we 
can assume without loss of generality that R g is continuous and has constant sign. If R g 
is negative, then J2*(l) < and hence \ g < 0. If R g > then clearly \ g > 0. But if 
Rg = 0,then J 2 * (1) = so X g = 0. If R g > 0, then 

Mf)> \\f\\ 2 Hl > ||/||^. 
and hence \ g > 0. On the other hand, if R g = 0, then J2* (1) = and hence \ g = 0. 

Finally, we show that Condition 2 implies Condition 3. If \ g < 0, then A(f,f) < for 
some smooth function /, so fi g < also. On the other hand, if \ g > 0, then A(f, f) > 
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for every smooth function / and /j, g > also. Since ||/||| 2 ~ I l/l 1^6 we conclude that 
J2* (/) ^ -hif) for every smooth non- vanishing function /. Hence if X g > 0, then fi g > 0. 
On the other hand, from previous implications we have shown that if ji g > 0, then \ g > 0. 
So if \ g = then /j, g — 0, which competes the proof. □ 

We conclude by showing that the sign of \ g can be detected if (M, g) has scalar curvature 
with R g > or R g < even when i? 5 is not continuous. 

Corollary 3.4. Suppose (M,g) is a compact Riemannian manifold with g e H s and 
s > n/2. If i? s < then X g < and X g = if and only if R g = 0. If i? 5 > tnen A 5 > 
and \ g = if and oniy if R g = 0. 

Proof: Suppose I? 5 < 0. Computing J2* (1) we find X g < 0. If X g = then (_R 5 , c) ^ M = 
for any constant c. Since R g < 0, it follows that for any smooth function / 

= (# 5 ,nrin/) > (R g J) (M s > (R g ,rimx f) =0 
and from Lemma 2.7 we conclude R g = 0. 

li R g >0 then it is clear that A ff > 0. If R g > and A ff = then fig = also. Let {/ n } 
be a sequence of non-negative functions in H 1 such that J2(fn) — > Aig = and such that 
1 1 /n 1 1 1,2 = 1. Then we have | |/ n | |^i is bounded and some subsequence converges weakly 
in H 1 to a function / such that 1 1/| \ L 2 = 1. From the weak lower semicontinuity of J2 we 
have -hif) = and hence V/ = almost everywhere. So / is a constant and moreover 
(R g , f 2 ) = 0. This implies (R g , c) = for any constant c, and using the fact that R g > a 
similar argument as before implies R g = 0. 

Finally, we note that if R g = 0, then Proposition 3.3 implies X g = 0. □ 

4. The Hamiltonian Constraint 

In this section we determine conditions under which the Lichnerowicz equation 

-aA g <p + R g <p = \a\ 2 g ^T 2 ^ - br 2 (f) K+1 (13) 

is solvable. Although we only require solutions when r is constant, it might be useful for 
future considerations of non-CMC rough solutions to treat here the case of a general function 
r. Hence we seek necessary and sufficient properties of the conformal data (M, g, a, r) 
leading to the existence of a positive solution of (13), including the case where r is not 
constant. 

A standard approach to finding a solution to (13) is to first make a conformal change to a 
metric with amenable scalar curvature and then to seek constant sub- and super solutions 
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0_ and + such that < 0_ < 0+. Having done this, a barrier argument proves that there 
exists a solution with < 0_ < < 0+. Propositions 6.2 and 6.4 of Section 6 show 
that this same barrier argument works for rough metrics. So we concentrate here of finding 
sub- and supersolutions. 

To illustrate the potential for difficulty in finding constant sub and supersolutions for rough 
conformal data, consider the case \ g > 0, and let us suppose for simplicity that the scalar 
curvature R g = 1. We then seek constants 0_ < + such that 

0+ > k^; 2k_3 - 1 

<P- < \(T\ 2 g(t>Z 2K ~ 3 -bT(f)1 +1 . 

If the conformal data are smooth, r is constant, and a is non-vanishing, then the existence of 
such constants follows easily from the asymptotic behaviour of the functions y i— > y- 2K ~ 3 
and y i — > — But if a vanishes somewhere we cannot find a positive constant such that 

0_ < laf^Z 2 ^ 3 -br^_ +1 

everywhere on M. Similarly, if a has singularities, then we typically cannot find a constant 
supersolution either. Since the difficulties posed by zeros are similar to those posed by 
singularities, we consider how the problem of zeros is treated for smooth conformal data. 

Suppose that (M, g) is smooth, R g = 1, r = and suppose a ^ has zeros. A solution 
to the Lichnerowicz equation in this case was found in [Is95] by finding a non-constant 
subsolution. Specifically, </>_ can be taken to be the solution of 

-aA 5 0_ = \cr\ 2 g A 



where A = max{l,maxM H ff }- An argument using the monotonicity of the function 



| 2-! 

naxji, maxjtf \<J 

y I— > y- 2K ~ 3 and the maximum principle then shows that 0_ is a subsolution. Although 
this argument cannot be adapted directly to the situation where a might have singularities, 
there is another way to understand the role that this subsolution plays. Suppose we ignore 
the constant A and simply solve 

2 

-aA g <pi + 0i = \a\ g . 

By the maximum principle we have 0i > 0, so we can use 0i as a conformal factor. Setting 
g\ = (f)\ K g we find that the scalar curvature of of g\ is 

Rgi = 0r 2 ^ 1 (-aA,0i + 0i) = 0^ \a\ 2 g = <t>T +3 N ffl 
where o\ = 0^ 2 o\ Solving the Lichnerowicz equation is then equivalent to solving 

-aA gi (f ) + R gi (l)=\a l \ gi 0~ 2k " 3 . 
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Since the zeros of <j\ are incorporated into R gi , finding constant sub- and supersolutions is 
easy. Moreover, there is no obstacle to finding constant sub- and supersolutions even if o\ 
has singularities. 

This motivates our approach for solving (13) when X g > 0. We start by making a conformal 
change to a metric with scalar curvature encoding both the zeros and the singularities of the 
conformal data, and then we find constant sub- and supersolutions. The specific choice of 
conformal change replaces specialized arguments used previously for subcases of X g > 
with a single unifying technique. 

Proposition 4.1. Suppose (M n ,g) is a compact Riemannian manifold with g e H s , 
s > n/2 and X g > 0. Suppose also that r is a scalar function in H s ~ l and a is a transverse- 
traceless tensor in if s_1 . Then there exists a positive solution e H s of the Lichnerowicz 
equation if and only if one of the following conditions holds: 

1. o ^0 andr ^ 0. 

2. X g > 0, o ^ and r = 0. 

3. X g = 0, o = andr = 0. 

Proof: We wish to solve 

-aA 5 + i? 5 0= \<j\ 2 g <T 2 *- 3 -br V K+1 - (14) 

From Proposition 3.3 there exists a conformal factor 0i e H s such that g>i = cf)\ K g has 
continuous scalar curvature R gi with the same sign as X g . If A 5 = 0, a = and r = 0, 
then 0i is a solution of (14) and we are done. So we can restrict our attention to Cases 1 
and 2 where o ^ 0. Letting <ti = (p^cr, it follows that there exists a solution of (14) if and 
only if there is a solution of 

-a^ gi <P + R gi <P= \a 1 \ 2 gx <\>- 2 ^ -br 2 ^ 1 . (15) 

From Proposition 6.1, since [R gi + 6r 2 )02 is non-negative and not identically zero, there 
is a solution 02 G if s of 

-aA ffl 02 + (^ 3l + jjr 2 J 2 = |<ti|^ 

Since o\ ^ 0, we have moreover from the maximum principles Lemma 2.9 and Proposition 
2.10 that 02 > 0. Let #2 = 02 K fi'i an d o"2 = 0^ 2 o"i- As before, there is a solution of (14) if 
and only if there is a solution of 

-aA, 2 + R 92 <P = \a 2 \ 2 - br 2 (f> 2K+ \ (16) 
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We are now in a setting where we can compute constant sub- and supersolutions. Since 

I |2 j2k+3 r 2 /-2k 
= M ff2 2 -0^ 2 > 

a constant + is a supersolution of (16) if 

I |2 <2k+3 7 2 j-2k1 / i ]2 i-2«-3 7 2 j2k+1 

Pick + such that 0^ +4 > max M ^ 2K_3 and such that 0^ > max M ^ 2K . This is 
possible since 02 G i? s and is hence continuous. Then 

|_ |2 .2 K +3i ^ |_ |2 j-2k-3 
F2| S2 02 0+>|^2| 52 + 

-6r 2 2 2 > + > -br 2 ^ +1 

and therefore + is a supersolution. Similarly, pick a positive constant 0_ such that 
02k+4 < min M 2 " 2K_3 and < 2k < minjvf 2 ~ 2k . Then 0_ is a subsolution and it is 
easy to see that 0_ < 0+. So from Propositions 6.2 and 6.4 proved in Section 6 there exists 
a positive solution G -£P of (16). 

To prove the converse direction, we need to show that if there is a solution with \ g > 0, 
and if at least one of r = or a = 0, then either A 5 > and <x ^ or A 5 = and cr and 
r both vanish identically. This is a consequence of Proposition 3.3. If a positive solution 
of (13) exists, we can set g' = 4> 2K g and a' = 4>~ 2 a to obtain 

R 9 > = W\]' ~ br 2 - 

If a and therefore a' = 0, then Corollary 3.4 implies A ff / < 0. Since \ g > we have 
A 5 = \ g i = and r 2 = 0. On the other hand, if r 2 = and a ^ 0, then Corollary 3.4 
implies X g > 0. □ 

Existence in the case \ g < is more delicate. Assuming the scalar curvature of g is negative 
and continuous, we want to follow the approach of Proposition 4.1 and solve an equation 
such as 

-aA 5 + R g <f) + 6r 2 0i = \a\ 2 (17) 

to make a conformal change to a metric with scalar curvature encoding the zeros and 
singularities of r and a. But R g < and we are not assured that we can solve this equation. 
On the other hand, if r ^ we can solve 

-aA g (f) + br 0i = \a\ 

as a substitute. Assuming that a ^ 0, then 0i will be positive and we can set g\ = 2K g. 
The resulting scalar curvature 



R 9 i = 0r 2K ~ 1 + ^0i - &r2 0i) 
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is problematic. Since R g is strictly negative the term R g — br 2 does not vanish when t 2 
does. This prevents us from finding constant supersolutions. On the other hand, if R g were 
a multiple of — t 2 then there would be no such difficulty. We now show that the existence 
of a conformal change to a metric with scalar curvature proportional to — r 2 is the precise 
condition required to solve (13) in the negative Yamabe case. 

Proposition 4.2. Suppose (M n ,g) is a compact Riemannian manifold with g e H s , 
s > n/2 and \ g < 0. Suppose also that r is a function in H s ~ l and a is a transverse- 
traceless tensor in H^ 1 . Then there exists a positive solution e H s of the Lichnerowicz 
equation if and only if there is a positive function 0i e H s such that g\ = vp 2fi g has scalar 
curvature R gi = —br 2 . 

Proof: In the case a = there is nothing to prove we see can assume that a ^ 0. First 
suppose that there exists a conformal factor 0i such that g\ = 4>\ K g satisfies R gi = —br 2 . 
Since \ g < this implies in particular that r ^ 0. Solving the Lichnerowicz equation then 
reduces to solving 

-aA gi ^ - br 2 ^ = \ai\ 2 gi ^- 2k -^ - br 2 ^ 2K+1 
Let <p2 be the unique solution in H s of 

2 2 

-aA gi (p2 + 6r 02 = Wi\ gi . 

Since r 2 ^ and since a ^ 0, the solution exists and is positive. Letting g 2 = ^ 2 K gi we 
have 

R 92 = <P 2 2 ^\-2br 2 ( j )2 + \a l \ 2 gi ) 

= -2br 2 r^ + \^\ 2 g A K+ "- 

We can now solve 

-aA g2 ^ + \a 2 \ 2 g2 - 25r 2 2 2 > = \a 2 \ 2 g2 ^~ 2 ^ - br 2 ^ +1 . (18) 
by finding constant sub- and supersolutions. 

Let + satisfy 2 f K+4 > maxM 4> 2 2hi ^ an d 0+^ > max^ 20 2 2k . Then, arguing as in 
Proposition 4.1 we have 0+ is a supersolution. Let 0_ satisfy 2k+4 < min^ 2 _2k_3 and 
02k < min^ 202~ 2k . Then 0_ is a subsolution with 0_ < 0+. Hence by Proposition 6.2 
or Proposition 6.4 there exists a solution to (18). 

For the converse direction, let us first make a conformal change to a metric with continuous 
negative scalar curvature via a conformal factor 0i. Let g\ = 4>\ K g and for future reference 
set o\ = <p 1 a. We wish to find a solution of the equation 

-aA gi <f> + R gi <l> = -br 2 (j) 2K+1 . (19) 
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To do this we seek non-constant sub- and supersolutions. Since there exists a positive 
solution of the Lichnerowicz equation, there also exists a positive solution + of 

-aA ffl + + = |<ti| 9i <j> 2 + K+S - br 2 ^ +l > -br 2 ^ +1 . 

So 0+ is a supersolution of (19). To find a subsolution we consider the family of equations 

-OA gi (j) e - R gi <p e = ~R gi ~ £T 2 . 

Since R gi < there is a unique solution to this equation for each e. When e = the 
solution is identically 1. So taking e sufficiently small, noting that convergence in H s 
implies uniform convergence, we can ensure that <p e > \ everywhere. Fix such an e. We 
claim that 0_ = r\<p t if r] is chosen sufficiently small. Indeed, pick r\ such that r)4> e < <p + 
and such that t] 2k < | min^ 0~ 2k_1 . Then 

-aA ffl 0_ + Rg^ = r]{R gi {2(t> e - 1) - er 2 ) 

< —rjer 2 

< -6r 2 (# e ) 2K+1 

< -br 2 ^ +1 . 

Hence 0_ is a subsolution with < (f)- < (f)+ and we conclude there exists a positive 
solution of (19). □ 

The apparently new observation of Proposition 4.2 is that the study of the Lichnerowicz 
equation in the Yamabe negative case completely reduces to the study of the prescribed 
scalar curvature problem. In the smooth category there exists an integral condition [Ra95] 
we can use to detect if we can make a conformal change to a metric with scalar curvature 
— br 2 . Suppose g and r are smooth, \ g < 0, R g is a negative constant, and Vol fl (M) = 1 
(since we can always make a conformal change to a metric with negative constant scalar 
curvature and Vol ff (M) = I, there is no restriction on g here other than \ g < and 
smoothness). Then there exists a conformal change to a metric with scalar curvature — br 2 
if and only if 

. , iMQ^jj dV g 

mf r9 * > -R g (20) 

fe* ff 2 dV g 9 

where 

= {/ G H 1 : / > 0, / ^ 0, and f M r 2 f dV g = 0}. 

Since condition (20) is well defined even for rough metrics, it is a good candidate for 
a necessary and sufficient condition even in the rough category. It remains to be seen, 
however, if this is true. We prove here a weaker result that suffices for finding CMC 
solutions of the constraint equations. 
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Proposition 4.3. Suppose (M n ,g) is a compact Riemannian manifold with g E H s , 
s > n/2 and \ g < 0. Suppose also that r a function in H s ~ l such that t > c for some 
positive constant c. Then there exists a positive function e H s such that g = 4> 2K g has 
scalar curvature Rg = —br 2 . 

Proof: We undertake a series of conformal changes designed to build a metric with scalar 
curvature that is negative, bounded away from zero, and singular wherever r is. We 
first make a conformal change through a function 0i given by Proposition 3.3 such that 
91 = 4>\ K 9 has continuous negative scalar curvature R gi . Next, let v G H s be the solution 
of 

-a A gi v — R gi v + r 2 v = -R gi 
which exists and is positive since R gi < 0. Setting 02 = 1 + v we have 

-aA ffl 2 + R gi (t)2 = {2R gi - r 2 )v. 

Hence the scalar curvature of #2 = <$i9\ is 

R g2 = 4> 2 2K ~ X v(2R gx -T 2 ). 

Since R gi , v, and 02 are all continuous, and since R gi is negative, we conclude there exist 
constants c\ > c 2 > such that 

-ci(l + r 2 ) < R g2 < - C2 (l + r 2 ). (21) 

To prove the proposition it is enough to find a positive solution of 

-aA ff2 + R 92 <p = -bT 2 <f (22) 

which we do by finding constant sub- and supersolutions. Pick constants 0_ and + such 
that 

< 2 _« < % 
b 

and 

A2« ^ £1 
b 



( 



1 

1 + — -, 



minfr 2 ' 



From (21) it follows that 0_ is a subsolution and + is a supersolution such that 0_ < + . 
Hence there exists a positive solution <fi E H s of (22). □ 

The uniqueness of any solution of the Lichnerowicz equation follows, as usual, from the 
fact that the functions y 1— > y^ 2li ^ s and y 1— > — y 2K are monotone decreasing. 
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Proposition 4.4. Suppose (M n , g) is a Riemannian manifold with g e H s and s > n/2. 
Suppose also that r is a function and a is a transverse-traceless tensor in if s_1 . If 0i and 
02 are solutions of (13), then either 

1. 0i = 02, or 

2. (j = 0, t = 0, A 3 = 0, and 0i = c02 for some positive constant c. 

Proof: Let gi = 2k , o\ = 4>^ 2 cr, and = 02/01- Then solves 

-aA fll + lo-il^ - 6r 2 = 0- 2 -- 3 - 5t 2 2k+1 . 



Hence 0—1 satisfies 



aA 5l (0 - 1) + H 2 - 6r 2 l (0 - 1) = \ ffl \ 2 (0- 2 ^ 3 - 1) - 6t 2 (0 2k+1 - 1). (23) 



'Si 



.91 



Since (0 — 1)( + ) e if 1 , and since — aA ffl (0 — 1) e if s 2 C H 1 , we can compute using 
Lemma 2.8 

(-aA^-lM*-!)"),^, = «.(v(*-l),V(*- l)W) (Msi) 

= f^aW-y.vOA-i))^ dv 9 . 

Hence 

f a(V(0-l),V(0-l)) ffl dK ff = 

(24) 

J l^l^ (0^ 3 - 0)(0 - 1) + 6r 2 (0 - 2k+1 )(0 - 1) dV g 

If > 1 then (0- 2 ^- 3 - 0)(0 - 1) < and (0 - 2k+1 )(0 - 1) < 0. So the integral on 
the right hand side of (24) must be non-positive. Since the integral on the left-hand side is 
non-negative, we conclude 

f a(V(0-l),V(0-l)) ffl dV gi =0 (25) 
f Iml^ (0~ 2k - 3 - 0)(0 - 1) + 6r 2 (0 - 2 * +1 )(0 - 1) dV gi = 0. (26) 

A similar argument using (0 — 1)( _ ) as a test function shows 

JT <1 a (V(0 - 1), V(0 - l)) gi dV gi = (27) 

f (0- 2K - 3 - 0)(0 - 1) + 6r 2 (0 - 2 * +1 )(0 - 1) ^ = 0. (28) 

From the continuity of — 1 and equations (25) and (27) we conclude is constant and 
therefore 0i is a constant multiple of 02- Suppose that 0i ^ 02, so ^ 1. Then from 
equation (26) or (28) we obtain 



51-0 



2 9 

and therefore |cri| = 0, r = and A 5 = 0. □ 
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We can now establish in the rough setting the necessary and sufficient conditions under 
which the Lichnerowicz equation is solvable when r is constant. 

Theorem 4.5. Suppose (M n ,g) is a compact Riemannian manifold with g E H s and 
s > n/2. Suppose also that r is a constant and a is a transverse-traceless tensor in H s ~ l . 
Then there exists a positive solution e H s of the Lichnerowicz equation (13) if and only 
if one of the following conditions holds: 

1. o ^0 andr ^ 0. 

2. X g > 0, o ^ and r = 0. 

3. X g = 0, o = andr = 0. 

4. X g < 0, o = and r ^ 0. 

In Cases 1, 2, and 4 the solution is unique. In Case 3 the solution is unique up to 
multiplication by a positive constant. 

Proof: The proof of the proposition in the case X g > follows directly from Propositions 
4.1 and 4.4. The existence of a solution if X g < and r^0 follows from Propositions 4.2 
and 4.3; its uniqueness is a consequence of Proposition 4.4. So it only remains to show that 
if X g < and if there is a solution, then r ^0. But this is a consequence of Corollary 3.4, 
arguing as at the end of Proposition 4.1. □ 

Remark. It is easy to extend these results to so-called scaled matter sources. In the presence 
of scaled sources, the Lichnerowicz equation becomes 

-aA g <P + R g <P = \<j\] <p- 2 ^ + p0-^ K - 6r0 K+1 (29) 

where p G H s ~ 2 and p > 0. The previous theorems carry over to this case by replacing the 
condition a = (resp. a ^ 0) with a = and p = (resp. cr ^ or p ^ 0). For example, 
the initial step of Proposition 4.2 when X g > and a ^ or p ^ is to solve 

-aA g (j) + for = |cr| + p 



rather than 



9 1 2 

-aA g (f) + br = |cr| 



We have only used features of the term |cr| that are shared by the term p0 
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5. The Momentum Constraint 

The treatment of the momentum constraint for rough metrics follows easily from the a priori 
estimates of Corollary 2.6 and standard arguments. We have the following properties of the 
vector Laplacian Al = div L. 

Proposition 5.1. Suppose (M n , g) is a compact Riemannian manifold with g E H s and 
s > n/2. Then Al : H s — > H s ~ 2 is Fredholm of index 0. The kernel of Al consists of 
conformal Killing fields in H s and the equation 

a l x = y 

with Y E H s ~ 2 is solvable if and only if 

for every conformal Killing field Z E H S (M). 

Proof: The estimates of Corollary 2.6 imply that Al is semi-Fredholm. On the other hand, 
the vector Laplacian of a smooth metric is well known to have index 0. Approximating g 
with a sequence of smooth metrics, using the fact that the index of a semi-Fredholm map 
is locally constant, it follows that the index of Al is also 0. If AlX = 0, then integration 
by parts yields = — (AlX, X)^ m g ^ = (LX, LX)^ M g ^ so X is a conformal Killing field. 
Finally, the conditions under which the equation AlX = Y is solvable follow from the fact 
that Al is Fredholm and self-adjoint. □ 

The York decomposition of symmetric, trace-free, (0,2)-tensors is a standard consequence 
of Proposition 5.1. 

Theorem 5.2. Suppose (M n , g) is a Riemannian manifold of with g E H s and s > n/2. 
Let y be the space of symmetric trace- free (0,2) tensors in H s ~ l . Then 

y = im L © ker div 

where L acts on H s vector fields and div acts on y. 

Proof: Suppose a E y. If X is a conformal Killing field, then (diver, X)^ Mg ^ = 
— (a, LX)( Mg j = 0. So from Proposition 5.1 there exists a vector field W E H s such that 
AlVT = div a. Since a = LW + (a — LW), and since div(<7 — LW) = 0, we have proved 
that y = imL + ker div. To show that the sum is a direct sum, it suffices to show that 
im L fl ker div is trivial. If a E imL D ker div, then a = LW for some W E H s . Since 
div a = 0, we have A^W = and hence W is a conformal Killing field. So a = LW = 0. 
□ 
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6. Sub- and Supersolutions 

We now turn to the existence theorem for the semilinear problem 

-A g u = F(x,u) (30) 

where F has the form 

F(x,y) = 2 i F j (x)G j (y). 

3=1 

We start with the following properties of the linearized operator. 

Proposition 6.1. Suppose (M n , g) is a compact Riemannian manifold with g E H s and 
s > n/2. Suppose also that V E H 8 - 2 and let Jzf = -A g + V. Then Jgf : H a -> /7 a - 2 
for any er e [2 — s, s] is Fredholm with index and we have the estimate 

\\ u \\h< 7 ^ 1 1 I \h ct - 2 + IMI# 2 - s - (31) 

The kernel of Jzf is a finite dimensional subspace of H s . IfV>0 and V ^ 0, then the 
kernel of Jzf is trivial and we have the estimate 

\\u\\h° < \ \^?u\\ H a-2. (32) 

IffE H a - 2 , the equation 

admits a solution if and only if (f, v)/ M g ) = v f° r every v E ker Jf. Finally, if u E H° for 
some a E [2 — s, s], and if Jzfu E H T ~ 2 for some t E [a, s], then u E H T . 

Proof: First consider the Laplacian A g . The estimates of Corollary 2.6 imply that A g is 
semi-Fredholm acting on H a with a E (2 — s, s]. On the other hand, the Laplacian of a 
smooth metric is well known to have index 0. Approximating g with a sequence of smooth 
metrics, using the fact that the index of a semi-Fredholm map is locally constant, it follows 
that the index A g is also 0. With the help of Lemma 2. 1 it is easy to see that Jzf is a 
compact perturbation of —A g and we conclude Jz? acting on H a with a E (2 — s, s] is is 
also Fredholm with index 0. Since Jzf is self-adjoint, it follows that Jzf is in fact Fredholm 
acting on H a with a E [2 — s, s]. The estimate (31) is trivial if a = 2 — s and follows from 
Corollary 2.6 for a E (2 — s\. Since the index of Jzf is 0, and since Jzf is self adjoint, the 
dimension of the kernel of Jzf acting on H s is the same as the dimension of the kernel of 
Jzf acting on H 2 ~ s . So these kernels coincide, which shows that the kernel of Jzf acting on 
H a for any a E [2 — s, s] is contained in H s ; the kernel is finite dimensional since Jzf is 
Fredholm. If V > and V ^ 0, then the weak maximum principle Lemma 2.9 implies the 
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kernel of Jzf is trivial and estimate (32) follows since Jzf is invertible. The conditions under 
which the equation 

Sfu = f 

is solvable now follow from standard Fredholm theory. Finally, suppose u G H a with 
a G [2 — s, s] and suppose 

St?u = f 

with / G H T ~ 2 where a < r < s. We want to show that u G H T . To see this, we note that 
since (/,t))/ Mj j = (Jzf w, v) r M g \ = for every v G ker Jzf we can solve 

Sfu = f 

for some u G H T and therefore u — u G ker J? C H s . Hence u G i/ r . □ 
The existence theorem for the equation 

-A g u=j]F J (x)G J (u). (33) 
i=i 

is a technical refinement of the constructive approach used in many previous papers treating 
the constraint equations [CB04][CBIY00][Is95]. In [Ma04] the technique was extended 
to asymptotically Euclidean manifolds (M n ,g) of class H s with s > n/2. That proof 
divided into two cases: 3 -manifolds and higher dimensional manifolds. We follow the 
same approach here and start with a result for the hard case of 3-manifolds that treats a 
more general class of coefficient functions Fj than those considered in [Ma04]. 

Proposition 6.2. Suppose 

1 . (M 3 , (?) is a compact Riemannian manifold with g G H s and s > 3/2, 

2. u-, u+ G H s are a subsolution and a supersolution respectively of (30) such 
thatu^ < u + , 

3. each Fj G H S ~ 2 (M) and is non-negative, 

4. each Gj is smooth on I = [inf(u_), sup(u+)]. 

Then there exists a solution u G H S (M) of (30) such thatu_ <u< u + . 

Proof: We first assume 3/2 < s < 2. Let 



V(x) = l+Y J F j 

3=1 



J\ x ) 



min G'a 
I 3 
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so that V G H S ~ 2 (M) is nonnegative and not identically zero. Let Fy(x, y) = F(x, y) + 
V(x)y and so that Fy is non-decreasing in y. By this we mean that if u, v E H s satisfy 
u > v , then Fy(x, u) — Fy(x, v ) > as a distribution. Finally, let Ly = — A g + V so that 
: iP — > iJ s_2 is an isomorphism. 

We construct a sequence of functions as follows. Let uo = u + , and for % > 1 let U{ be the 
solution of 

LyUi+i = Fy(x,Ui). 

Now 

Ly(u\ - Uq) < Fy(x, Uq) - F(x, Uq) - V (x)uq = 

since uo is a super-solution. Moreover, 

L v (ui - U-) > F v {x, u+) - FV(x, U-) > 

since u + > U-. From Lemma 2.9 we conclude uq>u\>u-. Now suppose -uo > u\ > 
• • ■ > u>i > U-. Then 

Ly(ui+l ~ u i) = F v (x, Ui) - F(x, < 0. 

Hence Ui + \ < Ui. Also, 

Ly (Ui+i - U-) = Fy(x, Ui) - F(x, U-) > 0. 

So Ui > Ui + i > U-. We obtain by induction for the entire sequence u+ = uq > u\ > 
u 2 > . . . > u _ m 

We claim the sequence is bounded in H S (M). From Proposition 6.1 we can 

estimate 

\\ui+l\\Hs{M) Z \\ F v(x,Ui)\\ H e-2^ M y (34) 

Now Fy(x, Mi) is a sum of terms of the form 

F(x)G(ui) 

where F e H S ~ 2 (M). Fix s' £ (3/2, s). From Lemma 6.3 proved below we have 

IIF^G^II^ < {\\G(ui)\\ L ~ + HG'^IIl-II^II^'I- 

Since m_ < u- L < u + , we have uniform estimates for each of the terms ||G('Ui)||Lcx> and 
1 1 G'(ui) 1 1 loo. Hence 

(35) 
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Now from interpolation we know that for any e > there is a constant C (e) such that 

IKII# S ' < C(e)\\ui\\ L 2 + e\\ui\\ H s. (36) 

Again, since w_ < U{ < u+, we have uniform estimates for | \ui\ \ L 2. Combining (35) and 
(36), taking e sufficiently small, we obtain 

\\ui+l\\H s {M) < dKH^CM) + C 

for some constant C independent of i. Iterating this inequality we obtain a bound for all % 

\\ui\\n s < ll^+llip + 2C. 



Hence some subsequence of (and by monotonicity, the whole sequence) converges 

weakly in H S (M) to a limit u^. 



It remains to see is a solution of (30). Now converges strongly to in H s for any 
s' < s and also converges uniformly on compact sets. Hence for any £ C°°(M), 

{F v (x, Ui) - V(x)u i+ i) dV g -> Moo)0 



So 



f (Vmoo,V0) dV g = f F(x, Uoo )<f)dV g , 

JM JM 



which shows that Uqo is a weak solution. So —AgU^ = F(x,u 00 ) in the sense of 
distributions. 

To handle the case s > 2 we use a bootstrap. First suppose 4 > s > 2. From the above we 
have a solution u in iJ 2 . Since 2 > n/2 = 3/2 and since 2 > s — 2 £ [0, 2], we know from 
Lemma 2.2 that c(x)f(u) £ H s ~ 2 . Since— A g u £ iP~ 2 , Proposition 6.1 implies w £ if s . 
We now obtain the result for all s > 3/2 by iteration. □ 

The following technical lemma completes the proof of Proposition 6.2. 

Lemma 6.3. Suppose (M n ,g) is a Riemannian manifold with g £ H s and s > n/2. 
Suppose also that F is a smooth function, a £ [—1,1], and s' £ (n/2, s). Then for every 
u £ H s and v £ H a we have 

\\vf{u)\\ H °{M) £ \\v\\h° [||/(«)IIl°°(M) + I \ f'i u ) I \l°°(M) \ \ u \ \h s '(M) ■ 
Proof: If a — 1, then 

\\vf(u)\\ H l {M ) £ \\vf(u)\\ L 2 {M] + \\V{vf(u))\\ L 2 {M) 

\vf(u)\\ L 2 {M] + \\Vvf(u)\\ L 2 {M] + \ \vf'{u)Vu\\ L 2 {M) 



< 

< \\r\ 



H\M)\\f(u)\\L°°(M) + \\vf'{u)Vu\\ L 2 M y 
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To estimate the last term we note that since v G H 1 , it is also in L p where 

1 _ 1 1 

p 2 n 

On the other hand, Vu e L g where 

1 _ 1 s' - 1 

q 2 n 



Since 



1 1 _ 1 1 1 s' - 
p q 2 n 2 n 



1 

= 2 + 

1 

<2 



'1 _^ 
2 n 



we conclude 



||u/'(u)Vu|| L 2 (M) < \\v\\ H i\\f'(u)\\ Lx{M) \\u\\ Hs/ . (37) 
This proves the result in the case s — 1. 
We prove the case s = — 1 by duality. If w e H l (M), then 

|(u/(«),w) (Mjfl) | < ||u|| ff -i||/(w)«;|| ff i (M) 

< 1 1«| [| |/(«) I U~(M) + I \f(u) \\ L °o(M) I |u| \ H s' {M) ] | kl |ffl(Af)- 

Hence 



\vf(u) 



\H- 



< 



|/H|| L oo(M) + ||/'H|| L oo(M)|| W ||^ (M) 



We have therefore obtained the result for s 
follows from interpolation. □ 



— 1, and the result for all s e [—1, 1] now 



The proof of Proposition 6.2 does not apply to manifolds M n with n > 3 since it depends, 
for example, on the fact that H 2 is an algebra. From the point of view of the constraint 
equations, this is not a serious drawback since the most interesting dimension is n = 3. 
Nevertheless, it would be nice to have an existence theorem in higher dimensions. This can 
be done following the approach of [Ma04]. The idea is to first obtain existence in the space 
W k,p where k is the largest integer with k < n and where | = \ — The full regularity 
in H s then follows from an interpolation argument and elliptic regularity. The proof from 
[Ma04] readily adapts to the compact setting using results and techniques from [CB04] in 
place of [Ma03] to obtain existence in W k,p . The bootstrap argument then follows exactly 
as in [Ma04] ; we omit the proof. 
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Proposition 6.4. Suppose 

1. (M n , g) is a Riemannian manifold with n > 3, g E H s and s > n/2, 

2. u-, u + E H s are a subsolution and a supersolution respectively of (30) such 
thatu- < u+, 

3. each Fj e H S ~ 2 (M) and is non-negative, 

4. each Gj is smooth on I = [inf(u_), sup(u+)]. 

Then there exists a solution u € H S (M) of (30) such thatu- <u< u+. 
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